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The gravitational interactions of the four-dimensional effective theory describing a general N- 
brane model in five dimensions without radion stabilization are analyzed. Both uncompactified and 
orbifolded models are considered. The parameter space is constrained by requiring that there be no 
ghost modes in the theory, and that the Eddington parameterized post-Newtonian parameter 7 be 
consistent with observations. We show that we must reside on the brane on which the warp factor is 
maximized. The resultant theory contains TV — 1 radion modes in a nonlinear sigma model, with the 
target space being a subset of hyperbolic space. Imposing observational constraints on the relative 
strengths of gravitational interactions of dark and visible matter shows that at least 99.8% of the 
dark matter must live on our brane in this model. 

PACS numbers: 04.50.-h, 11.25.Mj 



1. INTRODUCTION AND SUMMARY 



Over the past ten years, there has been a large amount of interest in extra-dimensional models of the Universe. 
Models such as the ADD model [TJ [5] and the Randall-Sundrum model [3J H] have demonstrated the possibility of a 
solution to the hierarchy problem through a modification of the fundamental gravitational scale, and the potential to 
provide interesting models of dark matter. 

Building on the success of the Randall-Sundrum model, many papers have considered various extensions to it, 
including bulk fields [5] , radion stabilization mechanisms [6] [7] , and models including more than one or two branes 
[8TfT5]. A wealth of knowledge of the phenomenology of these models has been accumulated [TUfT?] . A variety of 
techniques to analyze such models have been devised, such as linearized analyses [18H22] , the "covariant curvature 
method" [53], the "gradient expansion method" [2"M2"fj] . and derivative expansion methods However, relatively 
little work has gone into analyzing a general TV-brane model, and of the techniques mentioned previously, only 
linearized analyses are well equipped to analyze such a situation. Even then, linearized analyses require a background 
solution to perturb. 

In a previous paper [28] , we proposed an approximation scheme based upon a two-lengthscale expansion which can 
be used to evaluate a four-dimensional low energy action for five-dimensional braneworld models, and demonstrated 
its application to an uncompactified /V-brane model. In this paper, we analyze the physics of the model, starting from 
the four-dimensional action previously found, and also generalize the results to orbifolded models. Our motivation 
in analyzing general 7V-brane models is to determine whether the presence of extra branes may overcome some of 
the constraints the RS-I and RS-II models have, particularly with regards to radion stabilization requirements. We 
investigate the parameter space of the general model, and find regions in which the theory has no ghosts. The 
parameter space is further refined by imposing observational constraints from Solar System tests of gravity. We 
consider the possibility of placing dark matter and Standard Model fields on separate branes, and by comparing to 
observational data, find that the vast majority of the dark matter must reside on our brane in the models considered. 

This paper is organized as follows. We begin in Section[2]by describing the model and recalling the results from [2B] 
upon which we build in this paper. In Section[3] we diagonalize the field space metric of the radion modes for the entire 
parameter space. Only certain subsets of the model parameter space give rise to ghost-free four-dimensional theories; 
we derive the corresponding conditions in Section [4] In Section [5] we derive the physical four-dimensional low-energy 
action. Finally, in Section [6j we determine the observational consequences of this action, calculating the Eddington 
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FIG. 1: An illustration of the model a) before and b) after gauge fixing. The bulk cosmological constants A n , brane tensions 
a n , and metrics n g a p are labeled. 



parameterized post-Newtonian parameter 7 and the Newton's matrix for gravitational interactions between branes. 

The methodology considered here is also applicable to orbifolded models. In Appendix [XJ we show that the low 
energy theory for orbifolded models is very similar to that for uncompactified models. In Appendix [B] we discuss the 
spectrum of Kaluza-Klein modes in both orbifolded and uncompactified multibrane models. 



2. THE FOUR-DIMENSIONAL LOW ENERGY ACTION 



In this section, we recall the results from |28j that we use in the rest of this paper. 

Consider a five-dimensional model containing N four-dimensional branes, each with their own brane tension a n . 
We denote the n th brane by B n , where n is an integer with < n < N — 1. In between each brane there exists a bulk 
region of spacetime, which we denote TZq, . . . ,1Zn, with lZ n lying between branes n — 1 and n. In each bulk region 
lZ n we allow for a bulk cosmological constant A n . We allow matter fields represented by n (j> to reside on each brane, 
and include corresponding general four-dimensional matter terms in the action. The general five-dimensional action 
we begin with is 



r / f?(5)r 1 \ -^v— 1 » N—i 

S [ 5rs , V, n 4>] = / d 5 x^g J sJ - A(z r ) - E <W d 4 w n ^h+ £ n S m [ n h ab , >], (2.1) 

where grs is the bulk metric, T, S run from to 4, n x T = n x T (w^) gives the location of the n th brane B n in terms of 
the brane coordinates w% with < a < 3, k| is the five-dimensional Newton's constant, n h b is the induced metric 
on a brane, and A(x T ) takes the value A„ in the appropriate regions. 

Our first step is to partially fix the gauge so that the branes lie at y = 0, 1, . . . , N — 1, as seen in Fig. [l] 28 . 
Next, applying a separation of lengthscales technique, we find the general bulk metric which solves the high energy 
dynamics of the model, 

n ds 2 - e^ xC ^% b (x c )dx a dx b + X2y ^ V) dy 2 . (2.2) 
This metric is defined in each individual bulk region. The constants k n are defined by 



kn = i~ P' ( 2 - 3 ) 

and the function x( a;C : u) is continuous across all branes, although it has discontinuities in its first derivative in y. 
The solution requires that the brane tensions are tuned to some fixed values, 

where P n — sgn(x.j,) in lZ n (where P n is constant in each bulk region). 



3 



Having solved for the high energy dynamics, we now proceed to integrate them out of the theory, thereby generating 
a four-dimensional low energy action. We define the fixed parameters 



A n — 



1 



1 



1 1 



e„ = sgn 



(2.5) 
(2.6) 



for < n < N — 1. It is useful to note that e n = —sgn (a n P n P n +i) by Eq. (2.4). The values of the function x(a; c , n) 
evaluated on the branes become N scalar fields in the four-dimensional action, and we denote these by 



(2.7) 



where we use Xn — x( x0, i n )- There is residual parameterization freedom which implies that one of the fields \I/ n is 
nondynamical, but before we fix this freedom, we first give the four-dimensional low-energy action using the definitions 
so far. It is given by 



S[% b ^ n , n (, 



I 



d 4 x 



7 4^ 



N-l 



N-l 

E" 5 » 

n=0 



* 2 



lab, U 4> 



(2. 



where indices are raised and lowered using the four-dimensional metric i a b, and V a is the covariant derivative as- 
sociated with this same metric. In Appendix |Aj we show that an orbifolded iV-brane model gives rise to this same 
four-dimensional low-energy action with a rescaling of some parameters. Most of what follows from here is the same 
for orbifolded and uncompactified models. 

The residual parameterization freedom is 



X (x a ,y)^x(x a ,y)+S X (x a ) 

lab(x a ) -t labG 



-8x(x a 



(2.9) 
(2.10) 



under which the metric (2.2) is invariant. We can fix this freedom by specifying the value of x(x a ,n) for any n. In 
order to remain general, let us choose x(x a ,T) = 0, for some T with < T < N — 1. This causes the field ^>t to 
become non-dynamical. 

Some further field redefinitions now simplify the action. Let 



B n — 



A n 



(2.11) 
(2.12) 



Our dynamical scalar fields are now ip n , < n < N — l,n =/= T. Finally, we can define a four-dimensional effective 
Newton's constant as 



2^4 



1 

4k 



2 A T . 



(2.13) 



The action with these definitions is 



/ 4 / 





( N-l 


i? (4) [lab] 


i + E< 




\ n^T 




N-l 
n=0 


r n » n 
~B~labi 
_£>n 





N-l 



6 eTe«(V> n )(V a ^ 1 



(2.14) 



We next discuss the transformation to the Einstein conformal frame. Let P be the number of elements of the 
set {eT£n,0 < n < N — l,n=^ T} for which erf-n = +1, corresponding to the number of scalar fields with positive 
coefficients in the action. Note that < P < N — 1. Also, let M = N — 1 — P be the number of elements with 
£T£n negative, corresponding to the number of scalar fields with negative coefficients. It is convenient to relabel the 
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fields {ip n } based on which have positive kinetic coefficient (ipi, . . . , ipp) and which have negative kinetic coefficient 
(-0P+1, . . . , V'p+m), based on the action (2.14). We now define new coordinates Ci^ii ■ • ■ i@P-i an d r /i Ai, ■ ■ ■ , Aa/-i, 
such that 



(ipi,. • . ,tpp) = C (cos(0i),sin(0i)cos(0 2 ), ■ ■ ■ ,sin(0i) sin(0 2 ) ■ ■ ■ sin(0p_i)) , 
(V>p+i, • • • , iPp+m) = V (cos(Ai), sin(Ai) cos(A 2 ), . . . , sin(Ai) sin(A 2 ) • • • sin(A M -i)) . 



(2.15a) 
(2.15b) 



We choose r], £ > 0. All of the angular fields (0, and Xj) have a domain of to ir/2, as each \jj n is positive. We define 
the function 



N-l 



e = i + J2 ZT^i - 1 + c 2 - 



(2.16) 



n=Q 
n=£T 



and transform to the Einstein conformal frame using the conformal transformation g a b = 7 Q &|0|. 

Using these field definitions, the four-dimensional low energy action can be written in Einstein conformal frame as 



%afe, $ 



A n. 



d 4 Xy/^ge T sgn (9) 



R (i) [9ab] 1 



2k| 



7 AB ($ c ')5 ab V Q $ A V b $ i 



N-l 
n=0 



(2.17) 



Here, and V are associated with the metric g a b, $ A = (C> Vi • • • i 0p-i, Ai, . . . , Am-i), and 7ab( < I >C ') is the 
metric on field space, given by 



da = -f AB d§ d$ 



2 r 







-d( 







e 



(2.18) 



where e?S!p = cft^ + sin 2 (0i)<i02 + ... is the metric on the unit (P — l)-sphere, and similarly for dQ^. The parameter 
/i is defined by fi = %/6/«4. The coupling functions a„($ c ) are given by 



„2a„ 



1 

le|' 
lei B n 



0<n< N -l,n^T, 



(2.19a) 
(2.19b) 



where B n is given by Eq. (2.11), and ip n (<& c ) is defined by the relevant expression in Eq. (2.15) 



3. PARAMETERIZATION OF FIELD SPACE 



In this section, we find coordinates on field space which diagonalize the field space metric (2.18). We look at two 
special cases before analyzing the general case. 



A. Negative Definite Field Space Metric 



In the case M = 0, the general metric reduces to 



(i + C 2 ) 

M 2 



(3.1) 



This can be rewritten as 



da z = ~da z - ^ sin"* [ - ] dDr p , 



(3.2) 



where a = /xtan 1 (C), with < a < nfi/2. 
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B. Positive Definite Field Space Metric 



In the case of P = 0, the general metric reduces to 

.2 



1 



r, 2 dn 2 n 



1 - ?/ 2 

For the case where r\ < 1, this can be rewritten as 

dcr 2 = da 2 + ^i 2 sinh 2 



where a = /itanh with < a < oo. This is shown in Section [4] to be the only physically relevant case. 



For the case of 77 > 1, the metric (2.18) can be rewritten as 



da 2 = da 2 - fj 2 cosh 2 



dni 



(3.3) 



(3.4) 



(3.5) 



where a = /icoth (77), and < a < 00. 

We see that the two cases 77 > 1 and 77 < 1 are topologically disconnected, one being a metric on elliptic space and 
the other being a metric on de Sitter space, and so the divergence at 77 = 1 in the metric (3.3) is simply a coordinate 
singularity. 



General Case 



In the general case with M > 0,P > 0, the metric (2.18) is non-diagonal. It can be diagonalized using suitable 



coordinate transformations in the three different cases 0<O,O<0<1, and > 1 



/. 6 < 



For to be negative, we require from Eq. (2.16) that rj 2 — C, 2 > 1. Recall that 77 and £ are non-negative. We define 
new coordinates (a, b) by 



77 = a cosh 



C — a sinh 



where a > 1, b > 0. The metric (2.18) becomes 

,2 



da 1 



a 2 - 1 



db 2 



„, ^ Trda 2 + u 2 sinh 2 ( — 

or (or — 1) \n 

Defining c by a = cosec(c//j) with < c < 7T///2, the metric becomes 



d£li - v cosh z \~)d£E 



da 2 = sec 2 ( - ^db 2 + dc z + [i 2 sinh" ^ cosh z I - j dfi^ 



(3.6a) 
(3.6b) 

(3.7) 

(3.8) 



//. < e < 1 



In this regime, 77 > C, as previously, but with rj 2 — Q 2 < 1. We use the same coordinate definitions (3.6), but with 
< a < 1 and b > 0. The metric is the same as Eq. (3.7). This time, define c = /xsech _1 (a) with < c < 00, which 
gives 



da — cosech 



-db 1 -r «< :J /rsinlr ( dO 2 + /J 2 cosh 2 Q ) ,/<>; 



(3.9) 



as the metric. 
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in. i < e 



In this region of field space, £ > rj. We define coordinates (a, b) by 

b' 



rj = a sinh 
C = a cosh 



with domains of a > 0, b > 0. The metric (2.18) in these coordinates is 

,2 



a 

lTa 2 



-do + d6 - ^ cosh 



o?{l + a 2 ) 

If we define c = /icosech~ 1 (a) with < c < oo, the metric becomes 

b 



dfl 2 + n 2 sinh 2 



da 2 = sech 2 



dc 2 — fi 2 sinh 2 



cosh 



dnt 



dni 



(3.10) 
(3.11) 

(3.12) 
(3.13) 



The two cases < < 1 and > 1 are two coordinate patches on the same manifold. We see that the apparent 
divergence in the metric (2.18) at ?/ 2 — £ 2 — 1 is just a coordinate divergence; it delineates the boundary between 
topologically disconnected spaces (O > and 6 < 0). We show in Section [4] that only one of these cases is physically 
realistic. 



4. PHYSICALLY VIABLE MODELS 



In this section, we impose the constraint that all kinetic terms in the Einstein conformal frame have the correct 
signs, in order to exclude ghosts. This requires that the field spac e met ric h ave positive definite signature. Of the 
field space configurations, only those giving rise to the metrics (3.4 1 and (3.8) (with M = 1) meet this condition. We 



1 - P is the 
to have the 



investigate the constraints this imposes on the parameters of the model. 

Recall that P is the number of parameters in the set {ex£ n , n 7^ T} which are positive, and M = N — 
number which are negative. The metric (3.4) occurs when P = and M = N — 1. This requires all e r , 
same sign, except for ex which has the opposite sign. It also requires > 0. 

The metric (3.8) occurs with the correct signature when M = 1 and P = N — 2. This requires all e„ (including ey) 
to have the same sign except for one (not ex), which has the opposite sign. This metric also requires < 0. 

Combining these two cases, we see that all e„ (including ex) must have the same sign except one, which must be 
opposite. We now investigate what constraints these requirements impose. 

At brane B n , where the bulk regions n and n + 1 meet, there are four possible combinations for the parameters P n 
and P n +i, namely (P n ,P n+ i) = (—,—),( — ,+),(+,—) and (+,+). Furthermore, the bulk cosmological constant can 
either increase or decrease across the brane. The sign of the brane tension a n and the sign of e„ for each of these eight 
cases is given in Fig. [2j where the warp factor is plotted for each situation. Below, we refer to these eight possibilities 
as cases 1 through 8. We begin by looking at the situation where a single e„ is positive (0<n<iV— 1), and then 
look at the situation where a single e„ is negative. 



A. A single brane with e„ positive 



Recall that P n is the sign of the slope of the warp factor in lZ n . Using Pq = +1 and Pn = — 1 (which was assumed 
in deriving the four-dimensional low-energy action), we need a turning point in the warp factor somewhere in the 
progression of branes, which restricts us to either case 2 or case 6. Both of these cases have positive e, and so we 
require that all other e n are negative. Given that if the warp factor turns back upwards after turning downwards, it 
would need to turn around again using another case 2 or 6 which would introduce a second positive e, we see that the 
warp factor is only allowed to increase, turn around, and then decrease. The only way to continue increasing with 
negative e is using case 5, and the only way to decrease with negative e is using case 4. Thus, the progression of cases 
across the branes must go 



5,. ..,5,(2 or 6),4,...,4. 



(4.1) 
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< k n +i 


k n > &n+l 




1. 






5. 




(+,+) 






a -, e + 




a +, e - 
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FIG. 2: The behavior of the warp factor at a brane interface in the eight possible configurations. An increasing warp factor 
in a region has P„ — +1, while a decreasing warp factor has P„ = —1. In cases 2, 3, 6 and 7, the adjacent bulk cosmological 
constants can be equal. The horizontal axis in all plots is the y coordinate. 



It is unnecessary to have any branes with case 5 or 4 (the first or last case may be 2/6). Note that cases 2, 4, 5 and 
6 all correspond to positive tension branes. 

Given the growth and fall of the warp factor, there can only be one brane on which the warp factor is a maximum. 
We call this the "central" brane. Choose T to be this brane, such that x{ x °_jT) — 0, and so the warp factor is unity 
on the brane where the warp factor is a maximum. With the progression (4.1 ), e-r = +1, and all other e n — —1. We 



have P — and M = N — 1, and so we require that O > using these field definitions. 

We are interested in the sign of Q, to see if the requirement that > is met for the metric (3.4). As A n > 0, it 
is sufficient to know the sign of A^Q- We have 



A T Q = An 



n=£T 



(4.2) 



Now, given that the warp factor is a maximum on Bt and we know that P n = —1 for n > T, it follows that Xn > X«+i 
for n > T. Similarly, we have Xn < Xn+i for n < T. We now consider the expression for A n [Eq. (2.5)] based on 
what we know about P n and k n from the progression (4.1). 



A T = l/kT + l/k T +i, A n = l/k n -l/k n+1 (n>T), A n = l/k n+1 - l/k n (n<T) 
Thus, Q may be written as 



A T Q 



n<T n 



— e 



Xo 



E 

n>T 



k n +\ 



(e x " - e x " +1 ) + 



oXN-l 



N 



(4.3) 



(4.4) 



Each term in both sums is positive, and so 9 > 0. 

Thus, we see that a situation with all e n parameters negative bar one produces an action with no incorrectly signed 
kinetic terms. Furthermore, this choice of parameters requires all the brane tensions to be positive. Finally, the 
Ricci scalar in the action has positive coefficient, as er sgn(0) = +1. We investigate the properties of models in this 
parameter space in the remainder of this paper. 
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B. A single brane with e n negative 



Here, the number of possibilities is larger than in the previous case. By using the same logic as above, we find that 
the following progressions of cases are the only ways to meet the required conditions: 



Option 1 
Option 2 
Option 3 



l,...,l,5,l,...,l,(2or6),8,...,8 (4.5a) 
l,...,l,(2or6),8,...,8,4,8,...,8 (4.5b) 
1, . . . , 1, (2 or 6), 8, . . . , 8, (3 or 7), 1, ... , 1, (2 or 6), 8, . . . , 8 (4.5c) 



Each of these cases requires one or more negative tension branes. We consider each of these cases in turn. 
Option 1: 

Let the one negative e„ be ct, corresponding to case 5. One brane will have the maximum warp factor; call this brane 
X. Note that X ^ T, as brane T, being case 5, does not have the maximum warp factor. We now have ex — — 1, and 
all other e„ = +1, and so we have P = once again, which requires O > 0. Consider the sign of At@- We have 



A T G = A T -^A n e^. (4.6) 



We can once again calculate A n explicitly. 



A T = 1/fcr+i - 1/fcr, K = l/k n - l/k n+1 (0 < n < X - l,n ^ T), 

A x = l/k x + l/kx+i, A n = l/k n+1 -l/k n (n>X) (4.7) 
AtQ can then be expressed as 

X N-2 

A T Q = - — e Xo -Y^j- {e Xn - e^"- 1 ) - — e^- 1 - V (e x " - e*^ 1 ) . (4.8) 

n=l n N n=X n + 1 

Here, all bracketed terms are positive. Thus, O < 0, in contradiction of the requirement that 6 > necessary for this 
situation. 

Option 2: 

This case proceeds in exactly the same manner as Option 1, and we again find O < 0, in contradiction of the 
requirements for this situation. 
Option 3. 

This case is a little more complicated. Let T be the one brane with negative e, corresponding to case 3 or 7. Two 
branes will have a local maximum warp factor; let them be L and R (to the left and right of brane T) . Now, consider 
At&, which we require to be positive in this situation (as we once again have P — 0). 

A T e = A T - J2 A ne Xn - (4.9) 

This time, we have 

A n = ^- - — , < n < L, T <n< R, A n = — ^ , L <n <T, R < n, 

11 1 1 1 1 

k L k L+ i k T k T+1 k R k R+ i 

Combining these, we find 

A T e = - ^ - jr -L (e *» - e*»-) - >T J- (e*« - e^) 



A 1 ^ l e x«-i 

- I : jj- ^(e--^.,- (4.11) 

n=T+l n=R 

Once again, O is negative, and so this configuration also creates a contradiction. 
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FIG. 3: A diagram of the warp factor between branes and the associated bulk cosmological constants (dashed). Branes are 
represented as vertical lines. The bulk cosmological constants are negative, while the warp factor lies between and 1. 



C. The Effect of Negative Tension Branes 



From the above arguments, we see that the only ghost-free configurations are those which do not have any negative 
tension branes. This is consistent with the well-known local arguments for the instability of a negative tension brane. 
We note that by just using positive tension branes with the assumption that Pq = +1 and Pjy = — 1 (an d ignoring 
the requirement of the different e„ parameters having specific signs), the only possible combination is (4.1), and so it 
is the presence of negative tension branes which are giving rise to the instability. Any valid configuration which only 
has positive tension branes will not have this instability. 

The combination of cases (4.1) provides a rather tight restriction on the progressions of the bulk cosmological 
constant which can give rise to physically viable scenarios. Recalling that the bulk cosmological constants are negative, 
we require the bulk cosmological constants to increase across the branes monotonically to a maximum, and then 
decrease monotonically (see Fig. |3|. Note that in the special case where the first (last) brane has the maximum warp 
factor, then |A| can be monotonically increasing (decreasing). 



5. SPECIALIZING TO PHYSICALLY VIABLE CASES 



In this section, we specialize to the physically viable cases discussed above, and find a set of variables which simplifies 
the action. 



A. The Physical Action 



We previously found that the only physically viable configuration for the model is the configuration (4.1 ), in which 
the warp factor increases to a maximum, and then decreases again, with all brane tensions positive. We denote by 
n = T the index of the brane with the maximum warp factor, and call this brane the "central brane" . Specializing 
Eq. (2.14) to these parameters, we find 



S[% b ,ip n , n <t>] 



d x\J— 7 



24 



R (4) [%. 



N-l 

n=0 
n=£T 




N-l 



6]T(V a V>n)(V a VO 



n=0 



(5.1) 



This is the action in the Jordan conformal frame of the central brane. 
As P — 0, M = N — 1, the function Q is now given by 

N-l 



e = 1 - = 1 - 



V 



n=0 



(5.2) 
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and we know that > from the arguments of the previous section. We now follow the held redehnitions (2.15b) 
exactly, transforming into spherical polar coordinates. Let (Ai, . . . , Ajv_2) be angular coordinates such that 



V 
V 



cos(Ai) = f 
sin(Ai) cos(A 2 ) = /i 



(5.3a) 
(5.3b) 



1pT-l 

V 



sin(A!) . . . sin(A T _i) cos(A T ) = f T _ x 
sin(Ai) . . . sin(A T ) cos(A T +i) = fr+i 



(5.3c) 
(5.3d) 



4>N-2 
4>N-1 



sin(Ai) 
sin(Ai) 



. sin(AAr_ 3 )cOs(AAr_ 2 ) = fN-2 

. sin(AAr_ 3 )sin(A7v-2) = /jv-i- 



(5.3e) 
(5.3f) 



in the Einstein conformal frame, where g a b — ©7afc 

"i? (4) [s] 1 



Dehning a = /itanh 1 (77) with a > as in Section 3B we have our final four-dimensional low-energy action, written 
n the Einstein 

S[g a b, a, A„, n (f 



d xv — <? 



2«2 



, v N-2 (n-1 ~j 

- ( (V a a)(V a a) + M 2 sinh 2 ( - ) ^ \ J] sin 2 (A m ) }> (V a A„)(V a A„) 



n—1 K.m—1 



<"^ r I - ) 9ab, 4> 



N-1 



n=0 
n=£T 



sinh 



In a more convenient notation, the held space metric is 

da 2 = da 2 + ji z sinh^ 



B: 



dn 2 n , 



1 Tl 1 

y9ab, <P 



(5.4) 



(5.5) 



where dfi 2 = dX 2 + sin 2 (Ai)dA2 + ... is the metric on the unit N — 2 sphere. This is the metric on hyperbolic space. 

The target space will not be all of the quadrant of (N — l)-dimensional hyperbolic space for which all the held 
coordinates are positive, as we have yet to impose the constraint of having no branes intersecting, which was implicit 
in the derivation of the action. In the general case, these constraints are 



Xn < Xn+1, n<T, 
Xn > Xn+1, n>T, 



(5.6a) 
(5.6b) 



where \n is related to ip n by Eq. (2.12) 



B. The Effect of One Brane on Another 



Given the low energy action (5.4), it is interesting to ask about the effect one brane has on another, depending on 



how they are located. To investigate this, we consider two separate scenarios, one with iV branes, and one with N + 1 
branes, where an extra brane has been added after the last brane in the original scenario. The effect of this extra 
brane on rj 2 is to add an extra term to the sum (5.2 ). In the scenario with N + 1 branes, 



rf = V 2 +B N+1 e™, 



where 770 is the value of r\ in the scenario with N branes. 
The continuity of x(a: a ,y) across branes requires that 



dXn _ e XJV-i e ~2k N d N 



(5.7) 



(5. 
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where <1n is the geodesic distance between the now second last and last (newly added) branes. As exp(%jv-i) < 1 
(xt = is the maximum x)i this contribution to rj 2 becomes exponentially small as the distance to the new brane 
increases. Looking at Eqs. (I5.3I), we see that the change to the angular fields is also exponentially suppressed, and 



so the contribution of this new brane to the gravitational coupling is exponentially suppressed on all other branes. 
We therefore infer that the effect of the position of one brane on another, insofar as that information is coded into 
the radion fields, grows exponentially small as the distance between the branes increases. Given that the interbrane 
distances must be large compared to the AdS radii of curvature in order to meet the constraint from 7 (see Section 
6 A), this implies that the physics of a model with a large number of branes will dominated by the central brane and 



the nearest brane to it. 



6. OBSERVATIONAL CONSTRAINTS 



The theories (5.4) that are not ruled out by instabilities contain several massless radion fields, which will mediate 



long range forces and give rise to corrections to general relativity. Therefore, these theories will be subject to 
constraints arising from Solar System and other tests of general relativity. The nature of these constraints depends 
on which brane normal visible matter is assumed to reside. In this section, we investigate the extent to which these 
radion fields modify general relativity, and determine the corresponding observational constraints on the parameters 
of the theory. 



A. Eddington PPN Parameter 



The Eddington parameterized post-Newtonian (PPN) parameter 7, which measures deviations from general rela- 
tivity, is one of the most tightly constrained numbers from Solar System measurements of gravity. In this section, we 
compute this parameter from the action (5.4). 



As shown in Ref. [29] , for a theory of the form 



2 k 



N-l 



n=0 



' V 'gab, <f> 



(6.1) 



where § A are scalar fields and 7ab(3" C ) is the metric on field space, the Eddington PPN 7 parameter for observers 
on brane n is given by 



1-7 



where 



2 AB da n da 

n 

1 4 



n 2 _ „ . 

d$ A d<s> B 



and j AB is the inverse field space metric. For our theory (5.4), we have $> A = (a, Ai, 
is given by Eq. (5.5), and the functions a n are given by Eqs. (5.3) and (5.4). 



We calculate 7 for each of our branes. On the central brane, we find that 



t 2 



1 



(6.2) 
(6.3) 

. , Ajv-2), the field space metric 

(6.4) 



where r\ = tanh(a//i) has been used. As < 77 < 1, it is possible for T a\ to be sufficiently small on this brane to meet 
experimental constraints, which require that |30j 



|7- 1| < 2.3 x 10" 



(6.5) 



This constraint implies that the brane which is closest to the central brane must be at least 5 times the bulk curvature 
scale away from it (from Eqs. \2.12\ , fliT^ , and Eq. (5.18) in Ref. [28]). 

For the other branes, let 



p(n) 



n, n < T 
n — 1 , n > T 



(6.6) 
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in order to account for the hole in the sum over the matter actions in Eq. (5.4). For brane n, we calculate p a^ to find 



v 2 



> 



1 

3?7 2 

1 

3r/ 2 



1 + 




cot 2 (A,) 



=1 sin 2 (A r , 



+ (1 - 5 Pl N-2) 



tan 2 (A p+ i) 
n^ =1 sin 2 (A m ) 



(6.7) 



As < 77 < 1, none of these branes can give rise to a 7 parameter consistent with our observed Universe, and thus 
for this type of model not to be observationally excluded requires that we live on the central brane, where the warp 
factor is maximized. This implies that models of the form we are considering are unsuitable for explanations of the 
hierarchy problem, as no hierarchy can be obtained when considering Standard Model fields to be living on the central 
brane. Solving the hierarchy problem requires stabilizing at least some of the radion modes. 



B. Dark Matter Limits 



One of the motivations behind braneworld models is that the sequestering that occurs between matter on different 
branes may provide a natural explanation for the weakness of the coupling between normal matter and dark matter. 
Because of the different coupling factors of the metric to matter on different branes, there is a different Newton's 
constant for each brane, as well as different interaction strengths between matter on separate branes. As such, 
the Newton's constant becomes a Newton's matrix. In this section, we calculate the Newton's matrix measured by 
observers on different branes. 

The Newton's matrix depends on the brane on which the observer resides, since the units in terms of which the 
Newton's constant is measured vary from one brane to another. As the above section constrains normal matter to live 
on the central brane, we calculate the Newton's matrix from the perspective of the central brane. Generalizing the 



arguments presented in the appendix of |31j . for a theory of the form (6.1 ), we calculate the elements of the Newton's 
matrix to be 



M „2q, 



T cff 



8vr 



,AB 



da m da n 



(6.9) 



where G^p measures the strength of the gravitational interaction between matter on brane n and matter on brane 

„2 



Note that for m = n, the quantity in the brackets is 1 + n a§ 



When calculating the elements of (6.9), it is again convenient to write the quantities in terms of r\ — tanh(a//x). 
We also use p(n) (6.6), and similarly define q{m), in order to account for the missing term in the matter action sum 
in Eq. KM. We find 



TT 
eff 



G 



r'PP 



r<pi 



K 



i e 2a T 



Tp _ ^4 2a T 



87T 



^4 2a, 

8tt 



^4 2a, 
STr 6 




=1 llfc=l Sm ( A fe) 



+ (1 - Sp,N-2) 



tan 2 (A p+ i) 
nLiSin 2 (A, 



cot 2 (A,-) 



tnCiSin 2 (A fe ) ITLiSin 2 (A fe ) 



(6.10a) 
(6.10b) 

(6.10c) 
(6.10d) 



where m ^= n =/= T, and m < n. In all cases, the "1" in the outermost brackets arises from graviton exchange, while 
the remaining terms come from the exchange of scalar quanta. 

By considering the formation of the Sagittarius tidal streams, Kesden and Kamionkowski |32j have placed limits on 
the relative strengths of gravitational interaction between dark matter and normal matter. The constraint is roughly 



G 



M-DM 



\J Gm-mGdm-dm 



- 1 



5, 0.02 



(6.11) 
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where "M" indicates matter, and "DM" indicates dark matter. If we assume that all the dark matter lives on branes 
other than the central brane, we can calculate the constraints on our model that this provides, finding that 77 > 0.8. 
This disagrees with the constraint (6.5), which implies 77 < 6 x 10~ 3 . Thus, this model is unable to explain dark 



matter by positing the existence of matter fields on other branes 1 . 

We next consider the possibility that some fraction of the dark matter lives on our (central) brane, and some 
fraction lives on other branes. We can then calculate the percentage of dark matter which must reside on the central 
brane in order to be compatible with the observational constraints (6.5) and (6.11). On average, a mass M of dark 
matter will be composed of a mass aM on our brane, say, and (1 — a)M on other branes. The effective matter to 
dark matter coupling strengths will then be 



^eff 

/-iDD 
^eff 

^eff 



— G 

= Geff a 

= ^eff CX- 



eS 

TT 2 



-a) 2 +G^a(l-a) 
a). 



For simplicity, we use 



Gnn s-mvn 4 J2aj 

eff — ^eff ~ 7T~ e 
07T 



1 



(6.12a) 
(6.12b) 
(6.12c) 



(6.13) 



as the "off-brane to off-brane" coupling strength. Combining values for G^g,G^g and with Eqs. (6.12) in the 
constraint (6.11 1 and using rj 2 ~ 3.5 x 10~ 5 , we find a > 0.998, indicating that the vast majority of the dark matter 



must reside on our brane in this simplified model. 



7. SUMMARY AND CONCLUSIONS 



We have analyzed the observational constraints on a general five-dimensional braneworld model with arbitrary 
numbers of branes and without a radion stabilization mechanism, in the low energy, four-dimensional regime. The 
parameter space of such models was restricted by excluding ghost modes, and the phenomenology of the resulting 
models was analyzed. For such models to be viable, there is only one brane upon which Standard Model fields 
may reside, and such a configuration was unable to provide any benefit for the hierarchy problem, nor a natural 
explanation for the weakness of the coupling between normal matter and dark matter by sequestration. The Kaluza- 
Klein modes in such a model behave very similarly to the original RS-II model. Our model was not found to be ruled 
out experimentally, although observational constraints on the change in the value of GV between nucleosynthesis and 
today may do so. 

Overall, we found that models with N branes are quantitatively very similar to the two-brane case. Furthermore, 
uncompactified and orbifolded models were also found to be very similar, giving rise to identical four-dimensional 
low-energy theories. 



Acknowledgments 



We would like to thank S.-H. Henry Tye and Ira Wasserman for helpful discussions. This research was supported 
in part by NSF grants 0757735 and 0555216, and NASA grant NNX08AH27G. 



[1] N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, "The Hierarchy Problem and New Dimensions at a Millimeter," Physics 
\Letters B 429 (June, 1998) 263-272 hep-ph/9803315vl 



Note, however, that if the radion fields are stabilized, then it is possible to circumvent this restriction. As such, we can only rule out 
braneworld models with no moduli stabilization as an explanation for the observed weak interaction strength between dark matter and 
normal matter. 



14 



[2 

[3] 
[4] 
[5] 
[6] 

l<] 

[8] 

[9 
[10 

[11 
[12 

[13 

[14 

[15 
[16 
[17 
[18 

[19 

[20 

[21 

[22 

[23 

[24 

[25 
[26 
[27 

[28 

[29 

po; 

[31 
[32 
[33 



I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, "New dimensions at a millimeter to a fermi and 

superstrings at a TeV," \Physics Letters B 436 (1998) no. 3-4, 257-263] 

http://www. sciencedirect . com/science/article/B6TVN-3TVPMX6-5/2/8cda27091cc05185161631246f 4b6ef 6 



L. Randall and R. Sundrum, "Large Mass Hierarchy from a Small Extra Dimension," Physical Review Letters 83 (Oct., 
|1999) 3370-3373| [he^h/ 9905221 ~~" 
L. Randall and R. Sundrum, "An Alternative to Compactification," Physical Review Letters 83 (1999) no. 23, 4690-4693 
hep-th/9906064 



W. Goldberger and M. Wise, "Bulk fields in the Randall- Sundrum compactification scenario," Physical Review D 60 
|(Oct., 1999) 107505-9| 

W. Goldberger and M. Wise, "Modulus Stabilization with Bulk Fields," Physical Review Letters 83 (Dec, 1999) 
4922-4925, ~~ 

W. Goldberger and M. Wise, "Phenomenology of a stabilized modulus," \Physics Letters B 475 (2000) 275-279| 

F. Zen, Arianto, B. Gunara, and H. Zainuddin, "The Low Energy Effective Equations of Motion for Multibrane World 

Gravity," |hep-th/0511257v3| 

L. Cotta-Ramusino, Low energy effective theory for brane world models. PhD thesis, 2004. 
http : //www. gravity .port . ac .uk/theses/cotta-ramusino . pdf 

I. Kogan, S. Mouslopoulos, A. Papazoglou, G. Ross, and J. Santiago, "A three three-brane universe: new phenomenology 
|for the new millennium?," Nuclear Physics B 584 (Sept., 2000) 313-328||hep-ph/9912552| ~ 
I. Kogan, S. Mouslopoulos, A. Papazoglou, and G. Ross, "Multibrane worlds and modification of gravity at large 
scales," Nuclear Physics B 595 (Feb., 2001) 225-249, hep-th/0006030 



I. Kogan, S. Mouslopoulos, A. Papazoglou, and L. Pilo, "Radion in multibrane world," Nuclear Physics B 625 (2002) 



E. Flanagan, N. Jones, H. Stoica, H. Tye, and I. Wasserman, 


"A brane world perspective 


on the cosmological constant 


and the hierarchy problems," Physical Review D 64 (July, 2001) 045007, hep-th/0012129 




P. Brax and C. van De Bruck, "Cosmology and brane worlds 


a review,' 


Classical and Quantum Gravity 20 (May, 2003) 



R201 -R2321 

C. Csaki, "TASI Lectures on Extra Dimensions and Branes," hep-ph/0404096 
R. Maartens, "Brane- World Gravity," Living Reviews in Relativity (2004) 1-99. 
A. Perez-Lorenzana, "An Introduction to the Brane World," hep-ph/ 0406279 

J. Garriga and T. Tanaka, "Gravity in the randall-sundrum brane world ," Physical review letters 84 (Mar., 2000) 
|2778-81| [hep-th/9911055| 



H. Kudoh and T. Tanaka, "Second order perturbations in the Randall-Sundrum infinite brane-world model," Physical 
\Review D 64 (Sept., 2001) 84022| ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ 
H. Kudoh and T. Tanaka, "Second order perturbations in the radius stabilized Randall-Sundrum two branes 
\mode\," Physical Review L > 65 (MayT2002) 104034T - 



M. Carena, J. Lykken, and M. Park, "The Interval Approach to Braneworld Gravity," Physical Review D 72 (Oct., 2005) 
|084017[ |hep-ph/0506305vl| 



P. Callin and F. Ravndal, "Lagrangian formalism of gravity in the Randall-Sundrum model," Physical Review D 72 (Dec, 
2005) 10[ |hep-ph/0412109| 



T. Shiromizu, K.-i. Maeda, and M. Sasaki, "The Einstein equations on the 3-brane world," Physical Review D 62 (June, 
2000) 024012, 



S. Kanno and J. Soda, "Brane world effective action at low energies and AdS/CFT correspondence," Physical Review D 
|66 (Aug., 2002) 043526|~ 

S. Kanno and J. Soda, "Radion and holographic brane gravity ," Physical Review D 66 (Oct., 2002) 083506 

J. Soda and S. Kanno, "Holographic View of Cosmological Perturbations," Astrophys. Space Sci. 283 (2003) 639-644. 

T. Wiseman, "Strong brane gravity and the radion at low energies," Classical and Quantum Gravity 19 (Jan., 2002) 



|3083-3105[ |hep-th/0201 127| 

J. Bloomfield and E. Flanagan, "Four-dimensional description of five-dimensional N-brane models," Physical Review D 82 
| (Dec, 2010 J26] 



E. Flanagan, "The conformal frame freedom in theories of gravitation," Classical and Quantum Gravity 21 (Aug., 2004) 
|3817-3 829, gr-qc/0403063v3| 

C. M. Will, "The Confrontation between General Relativity and Experiment," Living Reviews in Relativity 9 (2006) 
no. 3, . 



R. Bean, E. Flanagan, and M. Trodden, "Adiabatic instability in coupled dark energy-dark matter models," Physical 
\Review D 78 (July, 2008) 023009, 0709.1128v2 



M. Kesden and M. Kamionkowski, "Galilean equivalence for galactic dark matter," Physical Review Letters 97 (Sept., 
|2006) 1313031 |astro-ph/0606566| 

T. Damour and I. K ogan, "Effective Lagrangians and universality classes of nonlinear bigravity, " Physical Review D 66 
| (Nov., 2002) 104024 



15 



Appendices 

Appendix A: Results on an Orbifold 

In this appendix, we derive the four-dimensional low-energy action of an orbifolded iV-brane model, and show that 
it is equivalent to the uncompactihed model up to the rescaling of parameters. 



Equivalent by periodic boundary conditions 



R-n+i , R-2 i R-i , R 1 R 2 

■aaH — —\ — 




B. N+1 B. 2 B.., I^Bq^J 

Equivalent by Orbifold 

FIG. 4: Diagram indicating how the branes are labeled in the construction of the orbifolded model. The orbifold symmetry 
identifies y with —y, and we impose the periodicity condition of identifying y with y + 2N — 2. To calculate the action, the 
model is broken up into 2(N — 1) bulk regions lZ n , but regions TZ n and lZ-„ coincide by the orbifold symmetry. 

We begin by describing the construction of the model. We use the notation established in Ref . [28 . Consider a 
model with N branes on an orbifold. The first and last branes are taken to be at the orbifold fixed points. The other 
N — 2 branes lie between these two branes on one half of the orbifold, and are duplicated on the other half by the 
symmetry. These regions lie on a circle, and so the coordinate describing the extra dimension will be periodic. To 
calculate the action for this model, we take there to be 2(N — 1) regions and 2(N — 1) branes. Let the first brane be 
labeled by Bq, situated at y = 0, where y is the coordinate describing the extra dimension. After gauge specializing, 
let there be N — 1 branes located at y = 1, 2, . . . , N — 1. In between the branes, we have N — 1 bulk regions. To 
account for the orbifolding, continue the extra dimension in the negative y direction, with another N — 1 branes 
located at y — —1, —2, . . . , — N + 1, with the coordinates y and — y identified. The y coordinate varies from — N + 1 
to N — 1, and these endpoints are identified under periodic boundary conditions in y. The branes labeled N — 1 and 
— N + 1 are thus the same brane. Labellings are described in Fig. [4] The action for this model is given by 




/np(5) \ N ~ 2 i r 

d 5 x ny /=^(— r -A n )+ £ - 2 \ d i w n V^h( n K + + n K-) 
]T aj ^^+X[V B ^j ]T WjX b /V] (A.l) 

n=-N+l \n=l n =-N+lJ 

The sums over branes which only run to N — 2 are written so because the branes — N + 1 and N — 1 are the same 
brane. Note that the brane tensions at the orbifold fixed points are included once only, while the brane tensions on 
the other branes are doubly included. This is just a choice of how to describe the brane tensions in the orbifold. The 
choice of the factor of 1/2 in the matter actions accounts for the doubling that occurs with the orbifolding. 

The procedure described in Ref. [55] may now be followed for each region. We gauge specialize to the straight 
gauge, before separating length-scales in the action. Writing the metric in each region as 

n ds 2 = e"*«^) »7 afc «, y n )dx a n dxl + "$ 2 «, y n )dy 2 n (A.2) 



with det(7) = —1, we can find the equations of motion at lowest order in the separation of length-scales. The 
following equations and boundary conditions arise, corresponding to Eqs. (3.12, 3.13, 5.3-5.7) in Ref. [35]. Note that 
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the equations in regions n and —n are identical, as required by the orbifolding condition: 



M<,») = " +1 xK,»), 
2 , n X 



Tt + 1 



n+1 



7a6 5?/ «) n )) 



= \ n i Clh n j bc ,y n i Cd n j d a,y - 3 ~ 2^"<f 2 A„, 

/ ™$ \ 

Tl 7 , = n 7 , ™7 6c "7 , - ™7 , 2 "y ^ , 

'oa.yt/ 'ao.jf / lcd,y lad.y \" A,y n ^ I ) 

1 "<I> 



= ^ n T b n i bc ,y n i Cd n i da ,y + ^y + "X 



y n. 



The boundary conditions at the first and last branes are 

— lab,y l yi= 0+ ' 



and 



= 



-Pi-k 5 ct = 



1 



AT-l 



7a6,y 



yjv _i=(iV-l)- 



1$ 



J/i=0+ 



iV-1, 



(A.3) 
(A.4) 

(A.5) 
(A.6) 

(A.7) 

(A.8) 

(A.9) 

(A.10) 
(A.ll) 

(A.12) 



(A.13) 



y N _ 1 =(N-l)- 



Equation (A.8 1 should be solved first. The solution (in matrix notation and suppressing indices n) is 

1 {x a , y) = A(x a ) exp (b(x°) j" $>(x a , y r )e- 2x( - x ^ y ">dy' 



(A.14) 



where A(a; a ) and B(x°) are arbitrary 4x4 matrices such that ■y has the properties of a metric. Combining this with 
Eqs. (A.5) and (A.6), we see that B is independent of region. The boundary conditions Eqs. (A. 10) and (A.ll) then 
imply that B = in all regions. Finally, the condition (A.5) then implies that A is independent of region, and so we 
can write n J a i,(x c , y) = "f a b(x c ) for all regions. 

The remaining equations of motion are then solved straightforwardly. Defining 

k n - 



we find 



and the brane tuning condition 



6 : 



\ y = 2P n k n n $ 



For the first and last branes, this condition is 



kiPi = - ^KgCTO, 



fcjV-i-Pjv-i = -Kfcrjv-i- 



(A.15) 
(A.16) 
(A.17) 

(A.18) 
(A.19) 
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The metric in each bulk region is 



n ds 2 = e n ^ xC ^% b {x c )dx a dx b + 



n x 2 y (x c ,y) 

4fc2 



dy\ 



(A.20) 



Following the prescription of [28], we now substitute this into the action ( A.l ) and integrate over the fifth dimension. 
The result is 



= / rf 4 



x\/ — 7 



AT-l 
n=0 



1 

2k§ 



'N-X 

E 

,n=l 



e Xn-l 



N-l 

y 

2 ^ 



oX„-l 



(Vx„-i) ; 



(A.21) 



where Xn(x a ) = n x(x a ,n). 

We now make the following definitions. 



A n — 


1 




1 


h P 
nj n ± n 


k n - 


hi-Pn+i 


A) = 


1 

k X P X 


1 

~~ kx 


A N -x = 


1 




1 


fejv-iPv-i 


k N - 



sgn 



1 



1 



eo 

e N-X 
#« 

With these definitions, the action is given by 

1 



d^x 



2«§ 



h P 

= sgn(-Pi) = 
= sgn(P A ,_ 1 ) 

= y/A n e^ 

'N-X \ 



kn+xPn+X 

-Px 
Pn-x 



N-X 



(A.22) 
(A.23) 
(A.24) 

(A.25) 

(A.26) 
(A.27) 

(A.28) 



^ (4) [Tab] E tn^l + 6 en(V a *„)(V a * n ; 



sn=0 



n=0 



N-X 

E 

n=0 



7/ C 



^2 

^7ab, > 



(A.29) 



This is identical to Eq. (2.8 1 above except for a factor of two multiplying l/4/c§, which arises from integrating each 
region twice rather than once. Otherwise, only the definitions of cq, Aq, £n-x and Aiv— l have changed, which corrects 
for the removal of the regions between the first and last branes and infinity in the bulk. Thus, the four-dimensional 



low energy action for this model is the same as for the uncompactified case (2.17), although some parameters have 
been modified. A spec ial case of the orbifolded model is the two-brane case, the Randall-Sundrum model [3]. In this 
case, the action (A.29) reduces to previously known four-dimensional actions [2"5] , 

Most of the analysis for the orbifolded scenario is identical to that for the orbifolded scenario. The only place 
time when the orbifolded scenario requires a separate analysis is when removing ghost modes. In the orbifolded case, 
we again want all e„ parameters to have the same sign except for one, which is opposite. Note that we now have 
6q = sgn(<7o) = —Pi and ejy-i = sgn(crjv-i) = P/v-i. For the first and last branes, we may only choose whether e is 
positive or negative, while for the intermediary branes, all of the previously discussed cases are possibilities. 

For a single positive e n , we need one of the following configurations: 



For a single negative e n , the options are 



-,5, 
-,4, 



.,5,(2 or 6), 4, 

•,4,- 

.,5,+. 

■ j 8, — , 

.,8,(3 or 7),1, 
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The analysis of each configuration proceeds exactly as in Sec. [4] We find that we must have a single positive e n , with 
all other e n negative. This implies that all branes must be positive tension, with the possible exception of the first and 
last branes, which may be negative. Again, the warp factor thus rises to a maximum and then falls again. If the first 
brane has the maximum warp factor, it has a positive tension, and similarly for the last brane. The four-dimensional 
low energy action specialized to such a configuration is described by (5.4) above. 

As the constraints on the Eddington 7 factor and the dark matter limits arise only from this action, the constraints 
on this orbifolded model are identical to those in the uncompactified model. 

In arriving at the four-dimensional low energy action (A. 29 1, we make the same approximations as in Ref. |28) . 
namely that the separation of length-scales is valid everywhere between the branes. However, we don't have any issues 
with the separation of length-scales breaking down towards infinity in the bulk, and nor do we need to invoke global 
hyperbolicity to constrain the behavior of the warp factor outside the collection of branes. Furthermore, the boundary 
conditions imposed by the orbifolding ensures that the degree of freedom B is projected out. In these regards, the 
orbifolded analysis is more robust than the uncompactified analysis. 



Appendix B: Kaluza-Klein Modes 



In this appendix, we venture away from the four-dimensional theory to investigate the Kaluza-Klein modes of our 
model. The methods and results here mimic the original RS-II model [4j closely. 

Consider an uncompactified model with N branes (with brane tensions tuned) and no matter. The solution for the 
five-dimensional metric will be 

ds 2 = e x(v) r] ab dx a dx b + dy 2 (B.l) 
where x y = 2fc n P n , and \ 1S continuous. Now consider metric fluctuations of the form 

ds 2 = (e^ Vab + h ab (x c , y)) dx a dx b + dy 2 . (B.2) 

Decomposing h ab into Fourier modes h ab (x c , y) = h ab %jj{y) exp(ip c x c ), where p c is a four-momentum with p 2 — —m 2 , 
we find to first order in h 



m e 



Our gauge choice is h a a = d a h ab = 0. Eq. (B.3) is equivalent to Eq. (8) in [4]. As discussed there, the solutions 
to this equation are Bessel functions (although here, they must be defined piecewise because of the piecewise nature 
of x)- There is a massless graviton mode, which has been integrated to give the four-dimensional effective graviton 



in our low-energy theory (5.4 1, and a continuum of massive Kaluza-Klein graviton modes, which in this paper were 
previously truncated. 

As in the RS-II model, there is no mass gap. Note that there are no so-called "ultra-light" [TOl El modes 
present in this model, as such modes occur in a model where the mass spectrum is quantized. Although the presence 
of extra branes complicates the mathematics, the physical effect of the Kaluza-Klein modes in our model is essentially 
the same as in the RS-II model. 

In an orbifolded model, the analysis of the Kaluza-Klein modes follows similarly, but the orbifolding condition 
implies that the mass spectrum is quantized, and we expect ultra-light modes to be present (see |33j and citations 
therein) . 



